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ABSTRACT 

We present a new method to quantify substructures in clusters of galaxies, based on the 
analysis of the intensity of structures. This analysis is done in a residual image that is the result 
of the subtraction of a surface brightness model, obtained by fitting a two-dimensional analytical 
model (/3-model or Sersic profile) with elliptical symmetry, from the X-ray image. Our method 
is applied to 34 clusters observed by the Chandra Space Telescope that are in the redshift range 
z G [0.02,0.2] and have a signal-to-noise ratio greater than 100. We present the calibration of 
the method and the relations between the substructure level with physical quantities, such as 
the mass. X-ray luminosity, temperature, and cluster redshift. We use our method to separate 
the clusters in two sub-samples of high and low substructure levels. We conclude, using Monte 
Carlo simulations, that the method recuperates very well the true amount of substructure for 
small angular core radii clusters (with respect to the whole image size) and good signal-to-noise 
observations. We find no evidence of correlation between the substructure level and physical 
properties of the clusters such as mass, gas temperature. X-ray luminosity and redshift. The 
scaling relations for the two sub-samples (high and low substructure level clusters) are different 
(they present an off-set, i.e., given a fixed mass or temperature, low substructure clusters tend 
to be more X-ray luminous), which is an important result for cosmological tests using the mass- 
luminosity relation to obtain the cluster mass function, since they rely on the assumption that 
clusters do not present different scahng relations according to their dynamical state. 

Subject headings: galaxies:clusters:general — large-scale structure of universe — X-ray:galaxy clusters 



1. Introduction 

Clusters of galaxies are the largest virialized ob- 
jects in the Universe, the upper limit of collapsed 
halo mass function. In an Universe dominated 
by a cosmological constant and cold dark matter 
(ACDM), dark matter halos are formed by gravi- 
tational instability from primordial quantum fluc- 
tuations in the mass density field. The amplitude 
of those fluctuations increases as they cease ex- 
panding with the Hubble flux, collapse and virial- 
ize, forming dense and relaxed structures. Smaller 
structures grow to larger ones through mergers, 
up to clusters of galaxies in the present time. In 
this hierarchical scenario of structures formation, 



clusters are thus dynamically young objects and 
contain e vidence of their recent past merging his- 
tory fe.g. iKauffmann fc Whit3ll993[ ). We can re- 



late substructures with the cluster dynamica l age 
fe.g. iRichstone et al.lflQoi : ISuwa et al.ll2003[ ): the 

more substructure (their total intensity) a cluster 
presents, the younger (dynamically speaking) it is. 

The hot intra-cluster plasma is a powerful X- 
ray source and its observation reveals the pro- 
jected spatial distribution of most of the bary- 
onic mass. X-ray studies of galaxy clusters are 
thus particularly relevant in this context, as they 
can give us clues to the dynami cal age of clus 
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((20051) - Abell 3921). Analysis of substructure 
in the intra-cluster plasma spatial distribution 
should help us determining the dynamical state 
of galaxy clusters. A very good review about the 
theory and observational status of the study of 
substructures based on X-ray data in clu sters of 
galaxies is given bv Ijeltema et al.l ( 20051 ). Here, 
we only briefly discuss some of the previous work 
on cluster substructures. 

Jones fc FormanI (1992) made the flrst X-ray 
systematic study of structures in galaxy clus- 
ters, visually analyzing 208 objects observed by 
the Einstein satellite, establishing that merg- 
ing must be a commo n phenomenon in clusters. 



Richstone et al.l (|1992f ) developed in an original 



theoretical study a relation between substructures 
and cosmology, where they put constraints on 
cosmological parameters by the fractional rate of 
major mergers in clusters. 

X-ray surface brightness allows us to per- 
form statistical tests such as cent roid and ellip- 



ticity variation ( Mohr et al. 19951 ). relating the 



dynamical ag e of cluster s with its morphology. 
iBuote fc Tsail (|l 9951 1 19961 ) developed a method to 
quantify X-ray substructures in clusters of galaxies 
from the moments of the expansion in Fourier se- 



ries of the X-ray surface brightness. iJeltema et al. 



(|2005[ ) used the same method, referred to as the 
power-ratio method^ in a sample of 40 clusters of 
galaxies observed by Chandra. They showed that 
clusters in general are less relaxed at z > 0.5, than 

Semi-analytic methods give an indication of 
the expected evolution of cluster substructure and 
its dependence on cosmological parameters, how- 
ever, the best method of constraining cosmological 
models is probably through the comparison with 
hydrodynamic clu ster simulations. For instance, 
Suwa et al. 1 (l2003h compared simulated clusters in 
a ACDM and an OCDM cosmology, at both z = 
and z = 0.5, using several methods for quantifying 
structure. They restrict themselves to comparing 
the ability of different statistical indicators in dis- 
tinguishing different simulated cosmologies, show- 
ing that cluster structure can potentially constrain 
l^A or the dark energy equation of state. 

Although a lot of effort has been done in order 
to advance our understanding about substructure 
in clusters of galaxies, from theoretical to numer- 
ical simulation studies, we propose in the present 



work a novel method of quantifying substructures 
that has a simple physical interpretation: the sub- 
structure level, the way it is deflned, reflects the 
fraction of the total X-ray luminosity that is emit- 
ted by the substructures, serving as a tool to un- 
derstand the underlying physical processes taking 
place during the cluster evolution. 

This paper is organized as follows. In section 2, 
we describe the sample selection and in section 3 
the data reduction and analysis are discussed. In 
section 4, the substructure level is deflned and our 
method is described, with its calibrations being 
discussed in section 5. In section 6 the results are 
presented and discussed and conclusions are flnally 
presented in section 7. The cosmology assumed in 
this paper is given by Qm = 0.3, Qa = 0.7 and 
Ho = 70 km s-^Mpc-^ 

2. Sample Definition 

Our method was applied to 34 clusters observed 
by the Chandra X-ray Telescope ACIS-I detector, 
with signal-to-noise ratios (S/N) greater than 100, 
and that are in the redshift range z G [0.02,0.2]. 
Figure [T] shows the cluster redshift distribution. 
There is an apparent gap within z G [0.12,0.14]. 
This is due to the incompleteness of the sample, 
but there is no particular redshift interval with an 
excess of objects and our results do not depend on 
the sample completeness. 

The criteria for clusters selection were chosen 
so as to ensure a suitable signal-to-noise ratio and 
a large enough image to work upon, without in- 
troducing bias for speciflc clusters. However, bi- 
ases that we do not control may affect our sample. 
Clusters are observed in time-competitive tele- 
scopes, so they must present something "special" , 
many times substructures and irregularities, that 
make them "worth" being observed. Therefore one 
should keep in mind this caveat, that it is possible 
that our sample may have a tendency to present 
more substructures than the average expected for 
all clusters in the redshift range z G [0.02,0.2]. 

3. Data Reduction and Analysis 

In order to obtain calibrated images without 
artifacts, adequate to be used with our method of 
substructure detection, it is necessary to follow a 
series of procedures of cleaning and flltering the 
X-ray data. Otherwise, we would have contam- 
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Fig. 1. — Histogram of the redshift distribution of 
clusters in our sample. 



ination that would be detected as spurious sub- 
structures. 

3.1. Data Reduction 

We have used the package CI AO 3.4. Initially a 
level 2 events file has been generated from a level 

1 events file, using the standard pipeline proce- 
dural! and the calibration files, CALDB 3.3.0. Pe- 
riods with high particle background (flares) were 
excluded using the Ic-clean script. At this point, 
a re-binned image with pixels corresponding to 
16 raw physical pixels (4x4, which roughly corre- 
sponds to 2" pixels) is created from the new level 

2 event file, in the energy band 0.3 to 7.0 keV. 
Then, we produce exposure maps and use them 
to obtain flat images from which the source points 
are removed by filling circles around each source 
with a random Poisson sampling with the same 
distribution as found in a circular region close to 
the source. Finally, we fit a 2D analytical surface 
brightness model. 

3.2. Surface Brightness 

The surface brightness profile is the projec- 
tion of the plasma emissivity along the line of 
sight. We will assume two radial analytical 
profiles for the surfa ce brightness: the /3- model 
( Cavaliere fc Fusco-Femiano 1976|) and the Sersic 
( 



Pislar et al.lll997l : IPemarco et al.ll2003[ ). 



In order to take into account the ellipticity of 
the plasma emission we use the following standard 



coordinates transformation: 

x' = {x — xo) COS — {y — yo) sin 0; 
y' ^{x- xq) sin 9^{y- yo) cos 6; 



(1-6)2' 



where (xq, ^o) is the X-ray emission center coordi- 
nates, 6 is the position angle, and e is the elliptic- 
ity. 

The /3-model may now be defined as follows: 



S(r) = So 



1+ I - 



-3/3+0.5 



b, 



where Tc is the core radius, P is the shape param- 
eter, and So is the central surface brightness. The 
parameter b corresponds to the background, and 
is supposed to be constant throughout the image 
(hence the importance of the exposure map cor- 
rection). 

The Sersic model is defined as follows: 



I](r) = So exp 



b, 



where a is the scale parameter, u (often repre- 
sented as 1/n) is the shape parameter and b is 
again the background surface brightness. 

Once we have the image correctly processed, 
we fit a 2D surface brightness model to it using a 
standard minimum squares method, and ob- 
tain the residual image, which is going to be the 
starting point for substructure quantification. 

We fitted the P and Sersic models for most of 
the clusters, and in the case where both models 
were fitted, we chose to use the one that gave the 
smaller substructure level (see below how the sub- 
structure level is defined and computed). In prac- 
tical terms, this is the same as choosing the fit with 
the smallest x^- The 2D surface brightness model 
fitted for each cluster is presented in Table O 

4. X-ray Substructures 

Previous studies on ICM substructure have 



been done, ei ther qualitatively (I Jones fc Forman 




Lagana et al 
19921 : 



-"^httpi/Zcxc. harvard. edu/ciao3.4/threads/createL2/ 



different techniques. There is, however, no method 
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that takes into account the ratio between the num- 
ber of counts on the residual and on the original 
images, which will be referred to as the resid- 
ual flux method. We describe here this method 
to quantify the substructure on the intra-cluster 
plasma emission. 

4.1. Substructure Level 

We start by defining a threshold for the residual 
image in order to identify the pixels which had a 
number of counts statistically significant above or 
below (positive and negative residues) the 2D sur- 
face brightness fitted model at the pixel position. 
The threshold in each pixel was defined as the 
square root of the number of counts of the model 
in the correspondent pixel, i.e., the expected vari- 
ance. 

Then we quantify the substructure level by 
computing the ratio between the total number of 
counts of the residual and original images (tak- 
ing the absolute value of the negative counts in 
the residual image and treating them exactly as 
the counts in the positive regions - after selecting 
those (in absolute value) above the threshold). By 
construction, the substructure level, 5, is defined 
as: 

s^'^ — , (1) 

i=i 

where CI is the number of counts of the i-ih. resid- 
ual image pixel and Cl is the number of counts of 
the z-th image pixel and n is the number of pixels 
of the image. 

We defined the substructure level this way be- 
cause it has a direct physical interpretation: it 
reflects the fraction of the total X-ray luminosity 
provided by substructures. 

The statistical uncertainties in the substructure 
level were computed using Monte Carlo simula- 
tions as described in the section [5^ 



5. Calibration of the Method 

5.1. General Case 

We may write equation ^ as: 

n 

Y,\Cl{t)-M[xt\ 

S{t)^^ , (2) 

^{Cj -h)xt 

i=l 

where M' is the model fitted to the image, which 
is decomposed into the cluster surface brigthness 
model and a constant background, i.e., M' = M + 
h. 

The number of counts of the i-th pixel, for a 
certain exposure time may be written as: 

Cl{t) = P[{h + S} + S^)xt\, (3) 

where 6^, S] and Sf are the expected number 
counts in the i-th pixel for an exposure time of 
t = 1, in an arbitrary time unit, from the back- 
ground, primary cluster and substructures, respec- 
tively. P{x) is the random Poisson deviate of 
the expected value x. In the limit when x ^ t^ 
P{x)^t±t^/\ 

By injecting equation (|3]) in equation (|2]) and 
taking into account that the sum of counts of 
the model is equal to the sum of counts of the 
main cluster plus substructures, i.e, Xir=i ~ 
ElLi(^'+^^'). we have: 

n 

E 1^" [{b^ + Si + Sf) xt]- (M, + b)xt\ 

S{t) = ^ , (4) 

J2{Sl + Sf)xt 

where b is the mean background level, i.e., b = 
1 b- 

5.2. Long Exposure Time Observation 

We now consider the limit of a very long expo- 
sure time. In this case, P[{bi -\- Sl -\- Sf) x t] 
{bi + Sl + Sf) X t, so the equation (j4|) takes the 
form: 

n 

Y^\(Sj + Sf)-M,-{b-h)\ 
limS{t) = ^ n ■ (5) 

E(^' + ^') 

i=l 
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We may write the model as: Mi = Sj + Di, 
where Di is the deviation on the i-th pixel due to 
Sf (the presence of substructures will change the 
model fitted in the i-th pixel by Di). Now the 
equation ([5]) takes the form: 

n 

J2\Sf-D,-{b-h)\ 
lim Sit) = , (6) 

which is different from the ideal case, 

n 

S' = ^ , (7) 

T.(sl + sf) 

in which the substructure level reflects exactly the 
fraction of counts provided by the substructures. 
However, using Monte Carlo simulations (which 
will be discussed in section [5^ to introduce sub- 
structure on model images of the clusters of the 
sample, one may correct this effect by introduc- 
ing a normalization factor in equation (|6]) for each 
cluster, allowing us to better estimate the true 
substructure level and quantify the systematic un- 
certainties involved in this method. 

5.3. Short Exposure Time Observation 

We consider now the limit when we have a very 
short exposure time. In this case the Poisson noise 
dominates over the expected value, P{t) t^/^, 
so the numerator of equation dH is dominated by 
noise, that is, it scales with t^/^, while the denom- 
inator scales with t, so S{t) ex t^/^/t = 

This property of the substructure level leads to 
the question: what is the minimum signal-to-noise 
ratio required for the method to be applied? In 
order to answer this question we must create im- 
ages of a cluster with substructure and vary the 
signal-to-noise to analyse how the substructure 
quantification varies. With this in mind we cre- 
ated images of different signal-to-noise of a cluster 
generated by a /3- model, with /3 = 2/3 and core 
radius equal to 20 pixels (in an image of 500 x 
500 pixels). We added a substructure 15 pixels 
away from the center of the main cluster. It has 
the same /3 and half the core radius and central 



surface brightness of the main cluster. Figure [2] 
shows how the measured substructure level var- 
ied with signal-to-noise ratio for different thresh- 
old levels (0.5, 0.75, 1, 1.25, 1.5, 1.75, 2 times the 
square root of the expected value for each pixel of 
the model). The residual images were smoothed 
using a Gaussian kernel of 1 pixel, which corre- 
sponds to 2" . The dashed horizontal line corre- 
sponds to the actual substruture level, i.e., the ac- 
tual number of counts that is provided by the sub- 
structure, and we see that the measured substruc- 
ture level converges asymptotically to this value. 
However, it converges differently according to the 
threshold used, hence the necessity to calibrate 
the method to obtain the factor for each cluster 
(according to the Gaussian smooth and threshold 
used) that will correct the value measured. The 
value we chose was {S/N)min — 100, although 
the method could also be applied to clusters with 
worse signal-to- noise ratio, which of course would 
increase the uncertainties. Figure [3] shows the 
signal-to-noise ratio distribution of the sample, 
which contains only clusters with signal-to-noise 
ratio greater than 100. We also want to stress that 
Figure [2] depends on the setup of the cluster and 
substructures and that it is purelly illustrative to 
show the different behavior of the measured sub- 
structure level as a function of the threshold used 
and the signal-to- noise ratio. 



5.4. Monte Carlo Simulations 

For each cluster, we generated an image of 
the main component using the best fit analyti- 
cal model. Then, we populated the images with 
substructures having random positions and inten- 
sities. A constant background was also added and 
to all image components were added a white noise 
following a Poisson distribution. For each cluster, 
200 realizations were done. 

Once the simulated images were generated, the 
same procedure used to real cluster images was ap- 
plied for all simulated images. We thus obtained 
a distribution for the substructure level that we 
compared to the actual level of substructure that 
was input into the simulated images, which we 
have control. 

The substructures added to the analytical im- 
ages had surface brightnesses described by a {3- 
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Fig. 2. — Measured Substructure vs. Signal-to-Noise 
Ratio. All of the curves correspond to a Gaussian 
smooth of 1 pixel kernel, and from the top to the bot- 
tom, the thresholds are: 0.5, 0.75, 1, 1.25, 1.5, 1.75, 
2 times the square root of the expected value for each 
pixel of the model. The dashed line corresponds to 
the actual substruture level, i.e., the actual number 
of counts that is provided by the substructures. Col- 
ors are used for better visual distinction between the 
curves. 
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Fig. 3. — Signal-to-noise ratio distribution of the clus- 
ters. 



model, with core radii and central surface bright- 
ness intensities that could vary between 25% to 
75% of the modeled cluster, the exact value being 
determined by a random variable. The number 
of substructures could also vary from (i.e., no 
substructure) to 3. 

In order to show that basically the quality of the 
substructure quantification depends on the size of 
the cluster compared to the whole image and the 
signal-to- noise ratio, we present Figure H] which 
shows how the corrected substructure level com- 
pares to the true values, for different cluster config- 



urations, in which different synthetic clusters were 
created, with fixed /3 = 2/3, core radius spanning 
from 20 to 80 pixels (whole image is 500 x 500 
pixels) and signal-to- noise ratios varying from 100 
to 700. We see in the left bottom plots that when 
substructures are close to the center of the clusters 
the method does not give a reliable result since the 
substructure is incorporated into the model when 
the surface brigfhtness fit is performed. Therefore, 
small angular core radii tend to give better re- 
sults since the amount of substructure which falls 
within the clustercentric distance is small. 

First, we made a linear fit of the measured sub- 
structure level against the true substruture level, 
i.e. for each cluster we had a relation: Sm = 
a-\-bxST^ where Sm^ St^ cl and b are the measured 
and true substructure levels, linear and angular 
coeficients, respectively. Once the fit was done, 
the corrected substructure level was computed by: 
Sc = {Sm — o)/h. In Figure [6] we see the mea- 
sured substructure level plotted against the true 
substructure level on the left panel and then the 
correction plotted on the center panel. 

The error bars were determined from the points 
distribution shown in Fig. [5l Starting with the 
cluster corrected substructure level, Sc^ we de- 
fined a symmetrical region Sc i ^Sc (horizon- 
tal dashed lines in Fig. [5j) where we have at 
least 18 data point^ in each side with respect to 
St — Sc = Then the asymmetrical error bars 
correspond to the range of 68% of the points in 
each side separately (red points in Fig. [5]). 

In Figure [6] we present the calibration for some 
(four) clusters of the sample, which were chosen 
because they represent different levels of substruc- 
ture (from Abell 907 with Sc = 0.062 to Abell 
2163 with Sc = 0.155) and different core radii 
and signal-to- noise. This figure shows both the 
measured against the true amount of substruc- 
ture (left panel) and the corrected substructure 
level (central panel) and the method we used for 
computing the uncertainties on the substructure 
quantification (right panel) . In Figure [71 we ilus- 
trate the Monte Carlo simulation with a very small 
sub-sample of the images created to calibrate the 



^Few points could determine the uncertainties erroneously, 
whereas too many would use points corresponding to very 
different corrected substructure level, so we decided (em- 
pirically) to use 18. 
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Fig. 4. — Illustration of the sensitivity of the method. For each plot we have the corrected (Sc) against the true 
substructure level (St) for each of the 200 simulated images. Each plot corresponds to different signal-to- noise (S/N) 
and core radii (Re) simulated clusters, with Rc given in pixels. The simulations were performed in 500 x 500 pixels 
images (see Figure [7|. For comparison, the dashed lines represent the 1:1 relation between the corrected subtructure 
level against the actual value. 



method for Abell 85. On the top left we see its X- 
ray image, as observed by the Chandra Space Tele- 
scope, and its simulated images containing ran- 
domly distributed substructures. 



6. Results and Discussion 

Now that we have measured and corrected the 
substructure level and estimated the error bars 
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Fig. 5. — Illustration on how the uncertainties on the 
substructure level are computed. 



within 68% confidencel level using Monte Carlo 
simulations, we may look for correlations between 
the substructure level, as we defined, and physical 
properties of the clusters. 

For correlations to be correctly assessed, it is 
necessary to well understand how variables are re- 
lated. Linear regression is a fundamental and fre- 
quently used tool in astronomy and it may seem 
surprising that such a statistical procedure, appar- 
ently simple, may be complicated and controver- 
sial (see , e.^., Isobe et al.lllQQOl : iFeigelson fc Babu 
|l992; Hogg et al. 2Q1Q[ for reviews). Briefly, when 
the scientific question clearly asks how one vari- 
able depends on the other, it is more appropri- 
ate to use OLS(F|X), Ordinary Least Square - 
the least square fit of the function Y(X), to quan- 
tify how the variables are correlated, with Y be- 
ing the dependent variable. However, when the 
scientific question does not clearly identifies the 
dependent variable, then it is recommended the 
used of OLS (Bisector) which is the bisector be- 
tween the OLS(r|X) and OLS(X|r) fits, the last 
case representing the fit inversion with respect to 
the variables. 



With respect to the size of the sample, when 
the size is small {N < 50, X the number of data 
points), resamphng methods such as Jackknife 
or B ootstrap should be used (jFeigelson fc Babu 
I992I ) to fit the data and estimate the uncertain- 



ties. 

The correlation streng th is estimated by the 
Pearson coefficient (see, iRodgers fc Nicewander 
19881), where its absolute value resides between 
and 1, with 1 meaning total correlation and 
none. The interpretation of the correlation 
strength depends on the context. A correlation 
of 0.9 may be very low if we are verifying a phys- 
ical law with high quality equipments, but may 
be seen as very high in social sciences for exam- 
ple, where there are many contributions of compli- 
cated variables. In Table [T] we give the two-tailed 
null hypothesis signifi cance for eac h Pear son corre- 
lation coefficient (See IPress et al.l ( 1992 ) for more 
information on how it is computed). 

In our case, we used the OLS(y|X) to fit rela- 
tions between the substructure level and physical 
parameters, whereas we used the OLS (Bisector) 
for the scaling relations, all fits performed using 
the Jackknife resampling method. 

6.1. Clusters Parameters 

In table [2] we give the corrected substructure 
level (6'c), M500, X-ray luminosity, tempera- 
ture, redshift, the 2D analytical surface brightness 
model fitted, core radius (Re) and signal-to- noise 
of the 34 clusters of the sample. 

We computed the substructure level and the 
2D analytical surface brightness model, while the 
other parameters were obtained from the litera- 
ture (iReiprich fc BQhringerll2002l: ISun et al.l[20ol 



Chen et al.l 120071 : iMaughan et al.l l201l[ ). M500 
and X-ray luminosity were corrected to a Hub- 
ble constant of 70 km s~^ Mpc~^, and then X- 
ray luminosity was extrapolated to the bolomet- 
ric band (0.01 - 100 keV) using K-correction and 
the XSPEC 12.0 MEKAL model (Mewe-Kaastra- 
Leidahl plasma emission code), since literature 
values were given for different Hubble constants 
and energy bands. Redshift s were obtained from 
NED (NASA/IPAC Extragalactic Database) with 
uncertanties varying from 10~^ to 10~^, therefore 
as they are extremly small compared to other un- 
certainties we use, they are not displayed in Table 

El 
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M500 - T 
Lx - M500 
Lx-T 


M500 = 0.312+--^ X ri-64±"-"^ 
Lx = 0.28ltn^^ X M5oJ-^4±° i^ 
Lx = 0.029l°:°J9 X T3-2o±o.29 


OLS (Bisector) 
OLS (Bisector) 
OLS (Bisector) 


0.94 
0.89 
0.91 


< 10-^ 

< 10-^ 

< 10-^ 




ocaie Kelations (Mign 


Substructure Level) 






M500 - T 
Lx - -M500 
Lx-T 


M500 = 0.378+^65 X Ti-4'^±"-iu 

T _ n /19/1+0-145 ^ Ti/T 2.03±0.15 
Lx - U.424_o ;l08 X ^^500 
Lx = 0.058tn?7 X T3-02±o.i8 


OLS (Bisector) 
OLS (Bisector) 
OLS (Bisector) 


0.92 
0.88 
0.90 


< 10-^ 

< 10-^ 

< 10-^ 




Scale Relations (Low 


Substructure Level) 






S-{l + z) 
S - Msoo 
S - T 
S-Lx 


S = 0.036+^^7 X (1 + ^)S-«±i-95 

C _ n n91+0 010 y /M- 0.62±0.20 

S = 0.016+° °J? X tO-74±o.30 
S = 0.040tn}^ X ^^o.i7±o.io 


OLS Y(X) 
OLS Y(X) 
OLS Y(X) 
OLS Y(X) 


0.35 
0.49 
0.36 
0.26 


0.042439 
0.003265 
0.036499 
0.137540 



Substructure Level vs. Physical Parameters 



Table 1: Table with best fits. Teraperature is given in keV, mass in lO^^M© and luminosity in lO^^erg/s. 



6.2. Substructure Level vs. Physical Pa- 
rameters 

Keeping in mind the different statistical ap- 
proaches, for all the correlations between the sub- 
structure level and physical parameters the OLS 
(F|X) was used, since the substructure level may 
depend on mass, temperature and luminosity but 
these quantities should not be dependent on sub- 
structure, the way it was defined. For the scaling 
relations we used the OLS (Bisector) since tem- 
perature, luminosity and mass have complicated 
relations connecting them. For instance, mass is 
one of the quantities which determines the cluster 
temperature, but temperature is used to compute 
the mass. Luminosity is the observed quantity 
(flux and redshift), although it depends on tem- 
perature. Therefore complicated relations exist 
between them, which made us use the OLS (Bi- 
sector). On the other hand, concerning the data 
size, since we have used 34 data points we chose 
the Jackknife resampling method to perform all 
the fits. The results of the fits are presented in 
Tabled 

6.2.1. Substructure Level vs. Redshift 

Figure [8] shows the substructure level as a func- 
tion of cluster redshift. We see that there is a 
dependence between the substructure level and 
redshift (5 oc (1 + ^)5.45±i.95)^ although within 
2.8(7 we find no evolution at all in the substruc- 
ture level. The dependence on redshift may be ex- 



plained by the fact that nearby clusters fill a larger 
detector area compared to more distance clusters, 
and as explained in the calibration Section (§[5j), 
they tend to have the substructure level understi- 
mated, since substructures that lie within small 
clustercentric distances are incorporated into the 
surface brightness fit and are hardly quantified. 
Furthermore the Pearson correlation coefficient of 
0.35 shows us a weak correlation, which translates 
as no significant structural evolution of the gas dis- 
tribution. We also note the strong scatter of the 
data points in the redshift range z G [0.02,0.2], 
showing that we find clusters in very different dy- 
namical states in this redshift interval, from those 
highly symmetrical to the very disturbed ones (see 
Figures [9l and [TQ|) . Such a scatter may be related 
to the young (dynamically speaking) age of mas- 
sive clusters. Abell 4038, which has the small- 
est substructure level, has been considered in all 
fits, however we present in Figures [HI [HI [12] and 
[T3I dashed lines representing the fits excluding it, 
since one could ask how much influence it has in 
determining the slopes of the curves. 

Compu ting: a temperature map based o n Chan- 
dra data, iMarkevitch fc VikhlininI ( 2001 ) clearly 
showed that Abell 2163 cluster (see Figure [TO]) 
is undergoing a major merger, which explains its 
high substructure level {SL = 0.155t^;^24)- 

6.2.2. Substructure Level vs. Temperature 

We see in Figure [TT] that hotter clusters show 
the tendency to present more substructures, how- 
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Table 2: Clusters parameters. Core radiu s (Rc) is given in pixel units (who le image is 500 x 500 pixels). Lit- 
eratur e refe rences (R) are given by : RB - iReiprich fc Bohringe rl (|2002h . S - ISun et all (|2004l ). C - IChen et al 
(|2007h . M - iMaughan et all (|2011f ). f The mass is given within l/3Ryir, Rvir = 1-35 Mpc. 
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Fig. 8. — Substructure level evolution. Solid and 
dashed lines correspond to the fit using all 34 clusters, 
and excluding Abell 4038 (the cluster with Sc < 0.004, 
on the bottom left of the plot), respectively. The fit 
with ah clusters is given by SL = 0.036to oo7 x (1 + 
z)5-45±i-95^ where we used the OLS Y(X). The Pear- 
son correlation coeficient obtained was 0.35. 

ever the Pearson coefficient of 0.36 shows a weak 
correlation. The slope of 0.74 ± 0.30 shows a pos- 
itive corrrelation within 2.5cr, however the weak 
correlation does not allow us to state any firm 
conclusion concerning the intensity of substrutures 
and gas temperature enhancements. 

6.2.3. Substructure Level vs. Luminosity 

We see in Figure [12] that substructure is basi- 
cally independent of X-ray luminosity, the Pear- 
son coefficient of 0.26 being the lowest between 
the substructure level and the physical parame- 
ters. Furthermore, the substructure level is com- 
patible with no dependency at all with the X-ray 
luminosity within only 1.7 a. 



Fig. 11. — Substructure Level vs. Temperature. For 
explanation on the fits, see caption on Figure [51 The 
solid line is given by SL = 0.016l°;°J^ x kT^-^'^^o.so ^ 
where the fits performed was the OLS Y(X) and the 
Pearson coefficient obtained was 0.36. 



1 ' 1 1 ' 1 ' 




1 






i 


j 




1 J 


: 1 




1 


1 i 

t 






i 1 - 
\ ■ 

1 



0.8 2 4 8 20 40 80 

L,^, 10"^"^ erg s"' 



6.2.4' Substructure Level vs. Mass 

Figure [13] shows the substructure level as a 
function of the cluster mass. We see that more 
massive clusters show the tendency to present 
more substructures. The relation between the 
substructure level and the mass has the form 



Fig. 12. — Substructure Level vs. Luminosity. For 
explanation on the fits, see caption on Figure [8] The 
solid line is given by SL = 0.040lg;gl^ x L/-^^^°-^°. 
The fits performed was the OLS Y(X) . The Pearson 
coefficient obtained was 0.26. 
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S (X 7\^0-62±o.20^ with a positive correlation within 
3.1(7, however, different from the temperature, lu- 
minosity and redshift correlations with the sub- 
structure level, the Pearson correlation coeficient 
is the largest (0.49) between them, presenting a 
strong (Null Hypothesis Significance= 0.003265 - 
there is only ^ 0.3% probability of not presenting 
correlation) relation between the amount of sub- 
structures a cluster presents and its mass. 




Fig. 13. — Substructure level dependency on mass. 
For explanation on the fits, see caption on Figure 
E The solid line is given by SL = 0.02ll°;°J° x 
^500 ^^^°"^°' where the fits performed was the OLS 
Y(X) and the obtained Pearson coefficient was 0.49. 



6.3. Scaling Relations 



Chen et al.l (|2007l ) constructed two cluster sam- 
ples based on the intra-cluster plasma central tem- 
perature, and they concluded that cooling-core 
clusters have different scaling relations compared 
to non-cooling-core clusters. Using the same idea, 
we created two different groups based on the sub- 
structure level. We computed the mean and me- 
dian substructure level of the sample (illustrated 
in Figure We choose the median, which is 

more robust regarding extreme data points, as a 
division line between high and low substructure 
level. Therefore, clusters were separated into two 
sub-groups according to their substructure level 



compared to the median value of the whole sam- 
ple. Numerically, a cluster was considered highly 
substructured if its substructure level were greater 
than S = 0.069755, and a low substructure level 
cluster otherwise. 



Median (Solid Line) = 0.069755 
Mean (Dashed Line) = 0.049724 




Fig. 14. — Substructure level distribution of clusters. 




Fig. 15. — Mass vs. Temperature. In red (triangles) 
and blue (circles) are the high and low substructure 
level clusters, respectively. The red (dashed) and blue 
(solid) lines are the best fit for the high and low sub- 
structure level clusters, respectively. The fit of the low 



substructure clusters is given by M51 



0.378 



+0.079 



^^i.49±o.io^ where the fit performed was the OLS (Bi- 
sector) and the Pearson coefficient obtained was 0.92. 
On the other hand the high substructure level clusters 
fit is given by M500 = 0.312t°;°^^ x kT^-^^^^om ^-^j^ ^ 
Pearson coefficient of 0.94. 
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Fig. 16. — Luminosity vs. Temperature. See Fig- 
ure [15] for information on the different points and 
lines. The fit of the low substructure clusters is 
given by Lx = 0m8tom7 x /cT^ o^^o is, where the 
fit performed was the OLS (Bisector) and Pearson 
coefficient obtained was 0.90. On the other hand 
the high substructure level clusters fit is given by 
Lx = 0.029lg;SJ^ X kT^-^o^o-^^ with a Pearson co- 
efficient of 0.91. 



Fig. 17. — Luminosity vs. Mass. See Figure [15] for 
information on the different points and lines. The fit 
of the low substructure clusters is given by Lx — 
0.424lg/i^^ X M5og o3±o.i5 ^ ^^^^^ performed 

was the OLS (Bisector) and Pearson coefficient ob- 
tained was 0.88. On the other hand the high substruc- 
ture level clusters fit is given by Lx = 0.28llo o55 x 
M.nn^^^^-^^ with a Pearson coefficient of 0.89 



The cluster segregation in high and low sub- 
structure level shows that hotter clusters are dis- 
tributed equaly between high and low substruc- 
ture level clusters as well as the more massive and 
luminous clusters do not reside in a preferential 
group (high and low substructure groups). 

The Pearson correlation strength coefficients 
for the scaling relations of these different groups 
are very similar, with the correlations for the high 
substructure group being slighly higher (roughly 
0.01 above). The slopes are very similar also, with 
the excepction of the the M — T relation being 
1.23<j compatible, while the L — T and L — M re- 
lations are 0.53cr and 0.40cr compatibles, respec- 
tively. However, differently from Figure [15] in 
which the slope is similar and the intercept off- 
set between the curves is small. Figures [T6] and [TTl 
show a clear off-set (2.07(7 and 2.23(7 for L — M and 
L — T respectively) between the curves for high 
and low-substructure level clusters. It appears 
that given a fixed cluster temperature or mass, 
the luminosity expected for low-substructure clus- 



ters tend to be higher. A possible explanation 
may be that low-substructure clusters, therefore 
more relaxed ones, have had enough time for the 
gas to accomodate into the gravitational poten- 
tial and become denser, which enhances th e X-ray 
luminosity. Figure 6 of I Chen et al.l (|2007l ) shows 
a very similar effect, where their cool-core clus- 
ters present higher X-ray luminosities, for a fixed 
temperature, compared to non-cool-core clusters, 
which would advocate in favor of low-substructure 
clusters, and therefore more relaxed ones being as- 
sociated with cool-core clusters. 

The diflFerent scale relations between low-high 
substructure level clusters suggests that substruc- 
tures are an important factor to bias scaling rela- 
tions, therefore it may affect the mass of clusters 
determination and thus the mass function that is 
used to test cosmological models, e.g. using the 
Press-Schechter Extended model etc. 
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7. Conclusions 

We have developed a new method to quantify 
X-ray substructures in clusters of galaxies based 
on the ratio between the number of counts in 
the residual and original X-ray images. We cal- 
ibrated the method and then applied it to 34 clus- 
ters of galaxies in order to obtain the substructure 
level dependence with physical parameters, such 
as mass, temperature. X-ray luminosity and red- 
shift. 

The calibration was done using Monte Carlo 
simulations, which showed that the method recu- 
perates very well the true amount of substructure 
for small angular core radii clusters (with respect 
to the whole image size) and good signal-to-noise 
ratio observations. 

The substructure level high scatter (spanning 
from less than 1 to 16%) in the redshift range 
z G [0.02,0.2] shows that clusters are found in all 
dynamical states in the local Universe: from those 
relaxed to completely disturbed. 

We have not found any strong evidence of cor- 
relation between the substructure level and physi- 
cal properties of the clusters, gas temperature, X- 
ray luminosity and redshift. However, there seems 
to be a correlation between the substructure level 
and the cluster mass, which is given by its Pearson 
correlation coeficient of 0.49. For our sample of 34 
clusters it represents a probability of not present- 
ing correlation of the order of 0.3%. 

The distinction between high and low substruc- 
ture level clusters has shown to be interesting, 
since different scaling relations were found with 
these two sub-samples (they present an off-set of 
~ 2(7 - given a fixed mass or temperature, low 
substructure clusters tend to be more X-ray lu- 
minous), which is an important result for cosmo- 
logical tests which use the cluster mass-luminosity 
relation to compute the mass function. 

A practical application of our method would 
be the identification of clusters of very low sub- 
structure level. Such relaxed clusters would be 
ideal laboratories for studies where the equilib- 
rium hypothesis is of paramount importance. As 
an example, the work of Bertolami et al. (2007), 
on the interaction between dark matter and dark 
energy, made use of the Layzer-Irvine equation, 
which must hold for a system in virial equilibrium 
when there is no interaction in the dark sector. 



They applied their method to Abell 586, based 
on the analysis of Cypriano et al. (2005), which 
suggests that Abell 586 is indeed a very relaxed 
cluster. 

Finally, it is important to say that the method 
itself is interesting since it concerns a new way to 
quantify substructures in clusters of galaxies, with 
a very simple physical interpretation: it reflects 
the fraction of the X-ray luminosity provided by 
substructures. 
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A907 




Fig. 6. — Illustration of the calibration applied to some clusters of the sample. Left: Measured Substructure, Sm, 
vs. the True Substructure Level, St, for all 200 Monte Carlo simulations. Center: Corrected Substructure, Sc, vs. 
True Substructure Level for all simulations. Right: Substructure corrected vs. (True Substructure Level — Corrected 
Substructure ) for the simulations. We only show examples for 4 clusters. See text for details on the corrections 
made. For comparison, the lines represent the 1:1 relation between the measured and corrected subtructure level 
against the actual value. 
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Fig. 7. — On the top left corner the X-ray image of Abell 85 is presented along with some of its simulated images 
containing arbitrary substructure. 




Fig. 9. — Abell 4038, a symmetric cluster, with few substructures. The substructures are basically due to the ACIS-I 
chip gaps, which represent only a very small fraction of the total counts (Sc — 0.004). 
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